Ising Model on Networks with an Arbitrary Distribution of Connections 
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We find the exact critical temperature T c of the nearest-neighbor ferromagnetic Ising model on an 
'equilibrium' random graph with an arbitrary degree distribution P(k). We observe an anomalous 
behavior of the magnetization, magnetic susceptibility and specific heat, when P(k) is fat-tailed, or, 
loosely speaking, when the fourth moment of the distribution diverges in infinite networks. When 
the second moment becomes divergent, T c approaches infinity, the phase transition is of infinite 
order, and size effect is anomalously strong. 
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The Ising model is one of the immortal themes in 
physics. It is a traditional starting point for the study 
of the effects of cooperative behavior. Networks with 
complex architecture display a spectrum of unique ef- 
fects K|-0|, and so they are an intriguing substrate. The 
simulation of the Ising model on a growing scale-free net- 
work has demonstrated that it is extremely far from 
that on regular lattices and on 'planar graphs' 

In this Letter, we report our exact results and results, 
which, we believe, are asymptotically exact, for the ther- 
modynamic properties of the Ising model on the basic 
construction for 'equilibrium' random networks. These 
networks are the undirected graphs, maximally random 
under the constraint that their degree distribution is a 
given one, P{k). Here, degree is the number of connec- 
tions of a vertex. Correlations between degrees of ver- 
tices in such graphs are absent. In graph theory, these 
networks are called 'labelled random graphs with a given 
degree sequence' or 'the configuration model' The 
earlier interest was mainly in the percolation proper- 
ties of complex networks and the spread of diseases on 
them Jl0|-p7t , and most of the analytical results were ob- 
tained just for this basic construction (however, see Refs. 
psj p(i[ , where the Berezinskii-Kosterlitz-Thouless perco- 
lation phase transition was studied in growing networks). 

Our results demonstrate the strong effect of the fat tail 
of the degree distribution on the phase transition in the 
Ising model. The most connected vertices induce strong 
ferromagnetic correlations in their close neighborhoods 
at very large temperatures, and so their role is very im- 
portant. Surprisingly, we observe very strong effects of 
these vertices, even when they, at first sight, must be 
insignificant, namely, when the first and the second mo- 
ments of the degree distribution are still finite but the 
fourth moment already diverges ((fc ) — > oo). 

It is convenient to use the power-law degree distribu- 
tion P(k) oc fc -7 for parametrization. Then, (fc 4 ) di- 
verges for 7 < 5, (fc 3 ) diverges for 7 < 4, and (k 2 ) is di- 
vergent for 7 < 3. When (fc 4 ) < 00, the phase transition 
is similar to that in the Ising model on high-dimension 
regular lattices. However, for 7 < 5 its nature is quite 
different. As 7 decreases and the role of the highly con- 



nected vertices turns to be more important, T c grows and 
the phase transition becomes more 'continuous'. Below 
7 = 4 it is of the higher order than the second in Ehren- 
fest's terminology. Furthermore, at 7 = 3, T c approaches 
00, and the order of the phase transition is infinite. 

When (k 2 ) diverges (7 < 3), the effect of the highly 
connected vertices is crucial. In the infinite network, 
long-range magnetic order is not destroyed by any tem- 
perature, and finite-size effect is very strong. Formally 
speaking, we develop a theory for the infinite networks, 
but it also allows us to describe finite-size effect. Our 
main results are presented in Tab. |. 

Intuitive arguments. — The large networks that are 
studied in this paper have a tree- like 'local structure'. 
When we start from a randomly chosen vertex, and add 
its first nearest neighbors, second, . . . , n-th with all their 
connections, the resulting sub-graph is almost surely a 
tree. Of course, in a finite size network, there is a bound- 
ary for n, above which loops appear in such a sub-graph. 

Cooperative properties, which we study in the present 
Letter, are determined by this 'local' environment, where 
vertices have a quite different structure of connections 
than in the entire net. Interactions are transmitted 
through edges, from vertex to vertex. Hence, the fol- 
lowing characteristic is crucial: the distribution of the 
number of connections of the nearest neighbor of a ver- 
tex. In the networks under consideration, it is kP(k)/(k). 
Then, the nearest neighbors of a vertex have the average 
number of connections (the average degree) (k 2 )/(k), its 
second nearest neighbors have the same average degree, 
and so on. Notice that this value is greater than the av- 
erage number of connections for the entire network (fc), 
and it is much greater than (fc) if (fc 2 ) is large. Therefore, 
we estimate the critical temperature of the Ising model 
on the network using the formula T c /J — 2/ \n\q/ (q — 2)] 
for the Ising model on a regular Cay ley tree 21 with the 
coordination number q = (fc 2 )/(fc). The result is 
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We will show below that this naive estimate is exact 
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7 > 5, (A; 4 ) < oo 
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TABLE I. Critical behavior of the magnetization M, the specific heat 5C, and the susceptibility x m the Ising model on 
networks with a degree distribution P(k) ~ fc -7 for various values of exponent 7. r = 1 — T/T c . The right column represents 
the exact critical temperature in the case (fc 2 ) < 00 and the dependence of T c on the total number TV of vertices in a network. 



General solution. — Consider the Ising model on an a 
network with the Hamiltonian: H = —JJ2(ij) Si^j — 
H^2 i Si, where the first sum is over all edges of the 
graph, the second one is over all vertices, J > and H are 
the energy of the ferromagnetic interaction and magnetic 
field, respectively. Hereafter, we set J = 1. It is known 
that the regular Cayley tree is solved exactly by using 
recurrence relations plfl . As networks under discussion 
have a local tree-like structure, we apply this method to 
the Ising model on such nets. Consider spin Sq on a ver- 
tex with fco adjacent spins Si^, i — l,2...fco- Due to 
the local tree-like structure this spin may be treated as 
a root of a tree. We introduce 



9i,i(So) = 5Z exp 

Si=±l 



J2SiS m + S S 1 , i + Hj2Si)/T 

Mm) I 



(2) 



where T is temperature. The indices I and m run only 
over spins that belong to sub-trees with the root spin 
#1,4, including S M . Let x M = 9i,i(—)/gi,i(+), then the 
magnetic moment M of the vertex is 



Af = (e 2 ^-n^M)/(e 2ff/T + II^ 



(3) 



i=i 



The parameters x\^ describe the effects of the nearest 
neighbors on the spin Sq. In turn, x\^ are expressed in 
terms of parameters x%,i = g2,i{—)/g2.i(+), I = l,2...fci,i, 
which describe effects of spins in the second shell on spins 
in the first shell, and so on. The following recurrence re- 
lation between x n j and x n+ u holds 




where we introduce the function 

e (-l+H)/T + e (l-H)/T x 



,(1+S)/T 



{-l-H)/T 



X 



(4) 
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If a spin SVi+i,z is on a dead end, then x n +u — 1. Note 
that at H > 0, x n 1 ^ 1, while x n \ ^ 1 for H < 0. For 



H > 0, it is convenient to introduce x n j — exp(— h ni i). 
Here, at a given n, h n ,i are positive and independent ran- 
dom parameters, which play the role of random effective 
fields acting on a spin in the n-th shell from neighboring 
spins in the n + 1-th shell. Then, Eqs. (J3J) and (^) take 
the form 



M 



e 2g / r -cxp(-E£ifti,0 
e 2ff/T +exp( _ E feo i?liji) ' 

' k n<j -l 
1=1 









h n ,j = - In < 
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(6) 
(7) 



At dead ends we have h n+ u = 0. At H = in the 
paramagnetic phase h n ,i — 0, while in the ordered phase 
h n ,i ^ 0. Eqs. (||) and (^) determine the magnetization 
M of a graph as a function of T and H . We emphasize 
that these equations are valid for any tree-like graph. 

While deriving the recurrence relations, we started 
from some spin Sq and then made the recurrence steps 
along sub-trees. While solving the recurrence relations, 
we start from distant spins, i.e., from large n, and descend 
along sub-trees to the spin Sq. The recurrence steps Eq. 
(^) converges exponentially quickly. Therefore, we can 
set for the starting spins h sa . In the limit n — > 00 the 
parameter is the fixed point of the recurrence steps. 
The thermodynamic behavior is determined by this fixed 
point, which is reached from the neighborhood of h = 0. 

Let us average the magnetization M over the en- 
semble of random graphs with a degree distribution 
P(k). Introducing the distribution function of h n i, 

^ n {h) = {S(h — hn t i), and its Laplace transform & n (s) = 
f °° dhe~ sh "$> n (h) we obtain from Eq. (|J) the average 
magnetic moment 



poo 2H/T _ -ft f 



ds 

27T?' 



,sh 
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Here P{k) is the probability that the vertex has k con- 
nections. Equation (Q) gives the recurrence relation be- 
tween ^n(s) and \l/„ + i(s): 



2 



ds' 



(9) 



with \l/ n+ i(s) = 1 at a dead end. Here P(k)k/(k) is 
the probability that the the neighbor in the n + 1— th 
shell of a vertex from the n— th shell has k connections. 
We start from distant spins with ^„ +1 (s) k, 1 and large 
ri and make recurrence steps toward smaller n until we 
reach the fixed point. In the limit n — > oo, the recur- 
rence procedure converges to ^(s). The fixed point ^f(s) 
is a solution of Eq. (||) in which ^ n (s) and ^ n +i(s) are 
replaced by *(s). Then Eq. |j) with *i(s) — > $(s) 
gives the exact expression for (M). At H = in the 
paramagnetic phase we have ^(s) = 1. 

Let us find the exact critical temperature, T c . Con- 
sider a starting function *f/ n (s) = exp(— s6), where S 
is small: < <5 <C 1, s<5 <C 1. After the first m 
recurrence steps, we obtain \E' n _ TO (s) = exp(— s5f m ), 
where / = (k(k - 1)) (fc) -1 tanh(l/T). For^ / < 1, 
the recurrence steps lead to the fixed point 9(s) = 1, 
which corresponds to the paramagnetic phase. / = 1 
at a certain temperature T c which satisfies the condition 
(k(k - 1)) (fc) _1 tanh(l/T c ) = 1, so that at T < T c we 
have / > 1 and the recurrence steps lead away from the 
point 4"(s) = 1. One sees that T c is given by Eq. ([!]). 
Thus, the estimate (|l|) is exact. This result is valid for 
any degree distribution with (k 2 ) < oo. In this range, 
the average number of the second nearest neighbors is 
z 2 = (fc 2 ) - (k) [0 (zi = (k)), so 



tanh — = — 

T c z 2 



(10) 



At zero temperature we find the solution ^(s) = t s + 
(1 — t s )S St o, where S St o — 1 at s — and S Si q = at s ^ 0, 
t s is the smallest root of the equation t a = J^k P(k)k 
t*- 1 . This solution gives (M) = 1 - ^P(fc)^, which 
is exactly the size of the giant connected component of 
the network. This result is evident. Indeed, only vertices 
in the giant connected component have non-zero sponta- 
neous moment, which is equal to 1 at T = 0. 

Ansatz. — To study the critical behavior of thermody- 
namic quantities, we develop an ansatz in the spirit of the 
'effective medium' approach. Notice that the right-hand 
sides of Eqs. (0) and (Q) depend only on the sum of the 
independent and equivalent random variables h n> j. So, 
let us use the following ansatz 



^ hn,i 



kh + 0(k 



l/2\ 



(11) 



where h = (h n ,i) is the average value of the 'effective field' 
acting on a spin. The larger k the better this approxima- 
tion. Therefore, the most 'dangerous' highly connected 



spins are taken into account in the best way. With this 
ansatz, 



(M)=5>(fc) 



,2H/T 



-kh 



-kh 



(12) 



e 2H/T + - 
k 

Applying the ansatz ([ll]) to Eq. (0) yields 

h= -(fe)" 1 ^P(fc)fclny[e-( fe - 1 )' 1 ] = G(h) . (13) 

h plays the role of the order parameter. At H = 0, h = 
above T c and is non-zero below T c . 

Thermodynamic quantities. — Let us describe the criti- 
cal behavior of the thermodynamic quantities of the Ising 
model on the infinite networks. For this, one must solve 



the equation of state (13) near T c at H = 0. 

The case (k ) < oo. The expansion of G(h) over 
small h has the form: G(h) = g\h + g^h 3 + . . .. For 
brevity, we do not present exact expressions for the co- 
efficients c/i. Substituting this expansion into Eq. (|l3| ) 
determines T c and the order parameter h « ar 1 / 2 as 
a function of the reduced temperature r = 1 — T/T c , 
where a = [12 (k(k - l)) 2 /((fc) (k(k - l) 3 ) Tj 1 / 2 . Note 
that the critical temperature T c that follows from our 
ansatz coincides with the exact result (0). At small 
h and H = Eq. (n2h gives the spontaneous mo- 
ment (M) w (k) ar 1 / 2 /^. The magnetic susceptibil- 
ity can be calculated from Eq. Jl2]), by differentiating 
over H and then takinginto account the dependence 
of h on H from Eq. (|l3|). We obtain x{H = 0) w 

At T > T c the suscep- 
{T/T c - l)- 1 but 



({ky/(2(k 2 ){k(k-2))))T~ 
tibility has the same behavior x(0) 
with the double prefactor. 

In our ansatz, the average internal energy per spin, 



(E) 



/N (the average is thermody- 



namic and over the ensemble of graphs) , at H = is 



(E) 



1 



(k) coth(2/T) 



2sinh(2/T) 



p -h I p -(k-l)h 

E^( fc ) fc i ,- kh ■ (14) 



1 



Substituting here h w ar 1 / 2 , we find that the specific 
heat C = d(E) jdT has a jump at T = T c : 



AC = (fc(fc-2))(fc 2 )a 2 /(8(fc)). 



(15) 



This jump disappears as (fc 4 ) approaches oo. When 
(fc ) diverges, we consider power-law degree distributions, 
P(fc)ocfc"^. 

The case 7 = 5. The divergence of (fc 4 ) leads to a 
logarithmic singularity of the function G(h): G(h) = 
gih + g 3 h 3 hx(l/h) + . . .. Solving Eq. (|lf) yields 

/t,Af ~T 1/2 /(lnT- 1 ) 1/a , 8C~l/\n T -\ x^^ 1 (16) 

[for any 7, we have SC(T > T c ) = 0]. Note that the 
critical behavior of x is not changed. 



3 



The case 3 < 7 < 5. The function G(h) has two lead- 
ing terms: G(h) f» gih + g^h 7 ^ 2 . Solving the equation of 
state (|l^) yields 



KM 



,1/(7-3) SC _ r (5- 7 )/(7-3) 



X 



(17) 



Notice that for 7 < 5 critical exponents of the magne- 
tization and specific heat differ from the standard ones. 
For 7 > 4, exponent a of the specific heat (SC ~ T~ a ) is 
above — 1 . Hence the second derivative of the free energy 
over T diverges at T c . Therefore, the phase transition 
is of the second order in Ehrenfest's classification. For 
7 < 4, we have a = —(5 — 7)/ (7 — 3) < — 1 and the 
transition turns to be of a higher order. The order of the 
transition tends to infinity as 7 — > 3. Nevertheless, for 
7 > 3, the susceptibility obeys the Curie law. 

The case 7 = 3. Here, (k 2 ) diverges. Formally speak- 
ing, this leads to the infinite critical temperature for the 
infinite networks [see Eq. (Q)]. In any finite network, 
(k 2X j < 00, and the critical temperature is finite, although 
it may be very high, T c = (k 2 ) / (k) (see below). We con- 
sider temperatures, which are much less than this crit- 
ical temperature, but where /» <C 1, so T > 1. Using, 
for brevity, the continuum approximation for the degree 
distribution, we obtain G(h) « ((h) h/2T) ln[2/( (k)h)]. 
Then, 



h^(2/(k))e- 2T /^, 
5C~T 2 e- 4T ^ k \ X 



M w e 
T -i 



-2T/{k) 



T r 



(fc)lniV 



at 7 = 3 . 



(18) 



r c ~ (7 u 2)2 (fc)jV(3-7)/(7-D for2<7<3. (20) 
(3-7)(7-l) 

The first expression can be compared with the simulation 
of the Ising model on the Barabasi- Albert growing net- 
work (7 = 3) with minimal degree fco = 5, so (k) = 10. 
From Eq. @ T c w 2.5 In N follows. The simulation § 
yields T c « 2.6 In N — 3. Recall that our results were 
obtained for the completely uncorrelated network, and 
correlations in growing networks are extremely strong. 

Discussion. — Our results may be compared with those 
for percolation on such networks ||l6|] and the disease 
spread within them [Q. Of course, the problems are 
distinct but one finds a great resemblance (the super- 
stability of long-range order when (k 2 ) — > 00, phase tran- 
sitions of higher order, etc.). All these anomalous fea- 
tures are determined by fat tails in the distributions of 
connections. Our final results were presented for uncor- 
related networks. However, the quantitative agreement 
with the simulation Q| of correlated nets shows that they 
are applicable in much more general situations. Further- 
more, our analytical results can be easily generalized. 
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Without the continuum approximation, we have, instead 
of (k) in the exponentials, a constant that is determined 
by the complete P{k). Equation ( [l8"| ) describes the be- 
havior of thermodynamic quantities at modest temper- 
atures in the situation, where the phase transition is of 
infinite order and at 'infinite temperature'. Notice the 
paramagnetic dependence x °t i/T. Note that the tem- 
perature dependence M cx exp[— IT j (k)] coincides with 
the result of the simulation Jjj] for the Barabasi- Albert 
model (7 = 3) [see Fig. 1 (a) from this paper]. 

The case 2 < 7 < 3. Again T c for large networks is 
very high. Using the expansion G(h) w gh~ t ~ 2 /T, we 
find, in the range 1 <C T <C T c , the following behavior 

h,M ~ T- 1/{3 -^\ SC ~ T- (7 - 1)/(3 - 7) , x ~ T- 1 (19) 

[compare with Eq. (|l8|)]. Notice that the susceptibility 
keeps the paramagnetic temperature dependence. 

Finite-size effect. — Equation ([!]) shows that T c di- 
verges when (k 2 ) — ► 00. In finite networks, (k 2 ) is finite 
because of the finite-size cutoff of the degree distribu- 
tion. In scale-free networks, it is estimated as k CMt ~ 
fcoiV 1 /^ 7-1 ), where N is the total number of vertices in a 
network, kg is a 'minimal degree' or the lower boundary 
of the power-law dependence, and (k) w k (j — l)/(j — 2). 
Then, actually, repeating estimations from Ref. M (see 
more detailed discussions in Refs. (l^jl^]), we obtain 
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